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HYBRID  COMPUTATION  OF  LEFT  VENTRICULAR  PERFORMANCE 
USING  NONINVASIVE  TECHNICS 


INTRODUCTION 

The  relationship  between  left  ventricular  function  and  externally 
recorded  cardiac  cycle  intervals  (7,  17)  has  been  appreciated  for  more 
than  50  years.  The  externally  recorded  signals  include  the  phonocardio- 
grarn,  electrocardiogram,  and  carotid  pulse  tracing.  An  explicit  definition 
of  the  relationship  between  physiologic  and  left  ventricular  systolic 
intervals  has  been  provided  for  both  animal  (1,  12,  13)  and  human  (3,  14) 
subjects.  Martin  et  al.  (8)  demonstrated  a clear  agreement  between 
simultaneously  recorded  internal  and  external  systolic  time  Intervals 
(STI).  A number  of  studies  with  human  subjects  and  patients,  as  reviewed 
by  Weissler  and  Garrard  (15),  have  additionally  related  STI  changes  to 
various  mechanical,  physiologic,  and  pharmacologic  interventions. 

The  application  of  STI  analysis  to  physiologic  investigations  and 
clinical  diagnosis  has  been  relatively  restricted  because  of  inherent 
limitations  in  STI  recording  and  interpretation.  Lack  of  specificity, 
for  instance,  is  illustrated  by  the  observation  that  STI  changes  resulting 
from  congestive  heart  failure  (10)  are  markedly  similar  to  those  accom- 
panying head-up  tilt  (16).  Additionally,  the  standard  methods  of  manually 
identifying  and  computing  the  intervals  are  very  time  consuming.  Standard 
heart  rate  corrections  in  accordance  with  regression  equations  for  group 
data  (15)  subsequently  require  even  further  computation,  and  the  universal 
validity  of  such  corrections  for  diverse  conditions  is  suspect.  The 
avoidance  of  beat-to-beat  analysis  in  these  standard  technics  may  signifi- 
cantly reduce  the  information  acquired,  as  evidenced  in  studies  (5,  9) 
demonstrating  identifiable  instantaneous  interval  changes  in  atrial  fibril- 
lation. Additionally  the  latter  study  showed  that  the  beat-to-beat  changes 
conveyed  information  pertaining  to  severity  of  the  mitral  lesion.  The 
need  for  automated  STI  analysis  has  subsequently  been  identified  as  a 
means  of  enhancing  the  reliability  and  accuracy  of  STI  analysis. 

Several  digital  computer  programs  (2,  4,  6,  11,  18)  have  been  written 
in  an  attempt  to  reduce  the  computation  time  associated  with  STI  analysis, 
but  each  of  these  programs  has  notable  disadvantages.  The  computerized 
analysis  is  restricted  to  averaged  data  segments  and  thus  fails  to  account 
for  beat-to-beat  variations  (references  2 and  18  are  exceptions).  Further- 
more, extremely  clear  signals  are  required  for  analysis,  thereby  prohi- 
biting application  to  stress  testing  and  other  potential  conditions  wherein 
signal-to-noise  ratios  are  reduced.  Finally  all  existing  programs  (except 
reference  2,  which  requires  a dedicated  minicomputer)  utilize  small 
Invest Igator-selected  quantities  of  cardiac  cycles,  and  this  precludes 
continuous  long-term  analysis. 
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This  paper  describes  a scheme  for  computerized  analysis  of  STI  that 
remedies  some  of  the  foregoing  deficiencies.  Externally  recorded  electro- 
cardiogram, phonocardiogram,  and  carotid  pulse  signals  are  input  to  a 
unilateral  hybrid  system  (analog  computer,  analog-digital  converter,  dig- 
ital computer)  to  provide  an  automated  printout  of  the  standard  left 
heart  intervals  on  a beat-by-beat  basis.  The  technic  provides  rapid 
analysis  of  large  numbers  of  cardiovascular  event  complexes  during 
periods  of  controlled  rest  or  stress.  The  automatic  nature  of  the 
pattern  recognition  process  effectively  removes  observer  bias  and  pro- 
duces meaningful  beat-to-beat  information.  This  methodology  thus 
enables  the  determination  of  both  transient  and  long-term  changes  in 
left  ventricular  function  that  result  from  pathologic  or  physiologic 
stress. 


METHODS 

Continuous  recordings  were  made  from  ten  human  male  subjects  for 
baseline,  centrifugal  acceleration  exposure,  treadmill  exercise  exposure, 
and  recovery  periods.  The  following  methodology  was  utilized  at  the  USAF 
School  of  Aerospace  Medicine  in  all  cases. 

Bipolar  ECG  leads  were  placed  on  the  subject.  Signals  were  pre- 
amplified proximal  to  the  subjects,  displayed  on  a standard  cathode-ray 
tube,  and  recorded  on  analog  magnetic  tape.  The  phonocardiogram  was  band- 
pass filtered  in  the  preamplifier  stage  with  lower  corner  frequency  of 
70  Hz.  The  upper  corner  frequency  was  set  in  accordance  with  character- 
istics of  the  individual  to  optimize  the  presentation  of  the  aortic 
component  of  the  second  heart  sound  (S2)  and  the  first  high-frequency 
component  of  the  first  heart  sound  (S^).  The  microphone  utilized  a 
standard  contact  crystal  which  was  acoustically  insulated  and  positioned 
on  the  chest  over  the  left  sternal  border.  The  carotid  pulse  was  detected 
by  an  air-coupled  diaphragm  with  the  associated  pressure  transducer 
isolated  from  the  subject.  The  carotid  pulse  diaphragm  was  strapped 
around  the  neck,  with  a minimum  exposure  area  of  the  diaphragm  positioned 
over  the  carotid  artery.  Voltage  output  was  preamplified  and  band-pass 
filtered  from  0.1  to  20  Hz.  The  phonocardiogram  and  carotid  pulse 
signals  were  again  filtered  with  the  same  band-pass  range  before  display 
to  develop  a rolloff  of  12  dB  per  octave. 

The  events  identified  by  electronic  processing  for  detection  of 
systolic  times  (ST)  are  illustrated  in  Figure  1 and  are  defined  as 
follows: 

1.  Rj  and  R2,  the  R wave  being  the  most  dependable  electrical 
event  occurring  during  the  cardiac  cycle:  the  onsets  of  the 
R waves  are  used  as  fiducial  points  to  bracket  the  cardiac 
cycle . 

2.  Q,  the  onset  of  electrical  activity  preceding  ventricular 
contraction:  the  beginning  of  the  Q wave  if  a Q wave  is 
present,  or  the  beginning  of  an  R wave  if  no  Q wave  is  detected. 
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PHONOCAROIOGRAM 


CAROTID  PULSE 
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ELECTROCARDIOGRAM 


Figure  1.  Simultaneous  tracing  of  the  electrocardiogram,  phonocardiogram, 
and  carotid  arterial  pulse  illustrating  the  events  necessary 
for  detection  of  systolic  times  (ST).  Sj  and  S2:  the  first 
and  second  heart  sounds.  B:  upstroke  initiation  for  each 
cycle  of  the  carotid  pulse.  E:  end  of  the  systolic  phase 
of  the  cardiac  cycle.  Rj  and  R2:  R waves  bracketing  the 
cardiac  cycle.  Q:  onset  of  electrical  activity  pAceding 
ventricular  contraction. 
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3.  Sj  and  S2,  the  first  and  second  heart  sounds  of  the  phono- 
cardiogram:  the  first  high  frequency  component  of  suffi- 
cient amplitude  for  each  of  these  heart  sounds. 

4.  B,  upstroke  initiation  for  each  pulse  cycle  of  the  carotid 
pulse:  the  intersection  of  a line  fitted  to  the  ascending 
limb  of  the  pulse  and  a horizontal  line  that  runs  through 
the  minimum  point  of  the  pulse. 

5.  E,  onset  of  incisura:  the  fastest  negative  moving  portion 
of  the  dicrotic  notch,  marking  the  end  of  the  systolic 
phase  for  the  carotid  pulse. 

All  but  one  of  the  above  points  are  found  by  analog  computer:  B, 
the  upstroke  initiation  of  the  carotid  pulse,  is  computed  in  digital 
mode  because  a line-fitting  technic  is  used;  and  this  cannot  be  per- 
formed in  the  analog  mode. 

With  ST  determined,  systolic  time  intervals  (ST1)  are  then  computed 
for  each  cardiac  cycle  and  are  identified  as  follows: 

1.  Total  electromechanical  systole,  Q-S2,  computed  as  the 
interval  between  Q-wave  onset  and  initiation  of  the  second 
heart  sound. 

2.  Left  ventricular  ejection  time,  LVET,  computed  as  the 
interval  between  carotid  upstroke  initiation  and  dicrotic 
qptch  occurrence. 


3.  Q-Sj,  computed  as  the  interval  between  Q-wave  onset  and 
first  sound  initiation. 

4.  si“S2»  computed  as  the  time  between  1st  and  2nd  heart 
sounds. 

5.  Preejection  Period,  PEP,  computed  by  subtracting  LVET 
from  Q-S2. 

6.  Isovolumic  Contraction  Time,  IVCT,  computed  by  subtracting 
LVET  from  S1-S2. 

The  unilateral  analog/digital  technic  utilized  to  locate  the  cardiac 
events  is  diagrammed  in  Figure  2.  Heart  sounds,  EKG,  and  the  carotid 
pulse  analog  signals  are  reproduced  from  magnetic  tape  (Ampex  2000A)  at 
double  the  recording  speed  (7.5  IPS)  to  reduce  analog  processing  time. 

The  analog  signals  are  then  low-pass  filtered  to  remove  tape  noise  and 
trunked  to  the  analog  computer  (EAI  680) . Heart  sound  signals  are 
scaled  and  squared  to  achieve  the  proper  amplitude  necessary  for  event 
detection  and  simultaneously  attenuate  nonusable  signals  (see  Fig.  3). 
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Figure  2.  Block  diagram  of  unilateral  analog/digital  technic. 
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Figure  3.  Block  diagram  of  680  analog  computer  program.  Sawtooth 
pulses  are  generated  marking  the  occurrence  of  systolic 
times  (ST). 


5 


1 


The  squared  signal  initiates  a monostable  timer  that  controls  an 
integrator  with  constant  voltage  input  which,  in  turn,  produces  a 
sawtooth-shaped  pulse.  The  beginning  of  the  sawtooth  coincides  with 
the  first  high-frequency  component  of  the  first  heart  sound,  and  the 
pulse  is  reset  to  enable  similar  detection  of  the  second  sound  (Fig.  4). 
Note  the  absence  of  an  S2  sawtooth  for  the  noisy  pulse  because  of  lack 
of  discrimination  between  S2  and  the  noise  episode.  The  monostable 
timer  provides  noise  immunity  because  any  noise  in  the  input  during  the 
timer's  on  period  will  not  affect  the  sawtooth.  The  sawtooth  allows 
the  transfer  of  an  ST  from  analog  to  digital  mode  at  a lower  sampling 
rate  than  would  ordinarily  be  required  to  achieve  the  desired  time 
resolution.  The  beginning  of  the  sawtooth  marks  the  occurrence  of  the 
ST  to  enable  calculation  of  the  time  of  the  ST.  A least-squared 
linear  fit  is  computed  for  six  points  on  the  upslope  of  the  sawtooth 
above  a selected  threshold.  The  intersection  of  this  line  and  the 
baseline  marks  the  ST.  This  sawtooth  method  is  used  to  transfer  all 
of  the  ST  found  in  the  analog  domain  to  the  digital  domain. 

The  EKG  is  differentiated  and  separated  into  sections.  The 
derivatives  of  the  Q wave  and  R wave  are  used  to  initiate  monostables; 
these  monostables  control  the  start  and  stop  times  of  sawtooth  waves 
that  mark  the  occurrence  of  the  Q and  R waves.  Analog  technics  could 
not  easily  be  implemented  to  locate  initiation  of  LVET;  hence  a digital 
procedure  was  devised  to  identify  the  initiation  of  carotid  pulse  up- 
stroke. Because  of  the  variation  of  carotid  pulse  from  individual  to 
individual,  the  end  of  LVET  was  difficult  to  locate.  In  actual  testing, 

I the  LVET  termination  could  be  detected  most  easily  by  using  the  nega- 

tive going  time  derivative.  This  testing  was  accomplished  over  a wide 
range  of  pulses.  Five  different  variations  in  wave  shapes  were  observed 
to  account  for  the  majority  of  cases.  In  Figure  5 (I),  the  pulse 
includes  only  one  readily  identified  notch;  therefore,  a single  saw- 
tooth is  generated  marking  the  end  of  LVET.  In  Figure  5 (II,  III,  IV), 
two  detectable  notches  appeared  on  the  pulse  with  the  second  marking 
the  end  of  LVET.  The  analog  computer  generates  separate  sawtooth 
pulses  on  different  channels  for  each  notch  with  the  digital  program 
automatically  selecting  the  second  of  two  notches.  In  one  other 
variation,  wherein  the  pulse  has  more  than  one  notch  but  the  first  is 
the  correct  one,  the  input  parameter  of  the  digital  program  provides 
for  the  selection  of  the  first  sawtooth  (Fig.  5 (V)). 

As  illustrated  in  Figure  2,  following  the  identification  of  ST  by 
analog  computer  and  marking  by  sawtooth  pulses,  all  analog  channels  are 
digitized  at  a frequency  of  125  samples/sec.  This  sampling  rate, 
required  for  accurate  definition  of  the  carotid  pulse  events,  is  used 
on  all  channels.  The  data  are  stored  on  1/2"  digital  magnetic  tape  and 
batch  processed  on  an  IBM  360/65  digital  computer.  Sawtooth  waveforms 
are  decoded  into  ST  on  a beat-by-beat  basis. 
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Figure  5.  Five  different  variations  in  carotid  pulse  wave  shape  that 
were  observed  to  account  for  the  majority  of  cases. 
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The  digital  program  initially  identifies  two  consecutive  R wave 
sawtooth  pulses  and  computes  indices  for  these  two  events.  The  R wave 
Indices  are  used  as  fiducial  points  to  bracket  the  events  in  a single 
cardiac  cycle.  A search  is  then  initiated  for  the  Q wave  sawtooth  at 
a maximum  distance  of  64  msec  preceding  the  first  R wave  index.  If  no 
Q wave  is  found,  the  R wave  is  accepted  as  the  beginning  of  electro- 
mechanical systole.  A character  (*)  is  then  printed  on  the  right 
margin  of  the  printout  for  this  beat  to  signify  that  no  Q wave  was 
found . 

Sj  and  S2  sawtooth  pulses  are  contained  on  the  same  digital  chan- 
nel, requiring  timing  consideration  for  correct  identification  of  the 
sound  pulses.  Sj  is  required  to  be  located  within  the  initial  25%  of 
the  time  for  R-R  interval  occurrence  to  be  counted  as  a first  heart 
sound.  A search  for  the  second  heart  sound  is  limited  to  the  remaining 
75%  of  the  R-R  interval. 

The  end  of  LVET  is  transferred  to  digital  form  by  sawtooth  pulses 
on  two  digital  channels,  depending  on  the  individual  variation  in  the 
notches  of  the  carotid  pulse  (Fig.  5).  With  most  cases  wherein  saw- 
tooth pulses  are  recorded  on  both  channels,  the  second  channel  sawtooth 
is  used  to  mark  the  true  end  of  LVET  (Fig.  5 (II,  III,  IV)).  The  pro- 
gram automatically  selects  the  second  of  the  two  notches.  In  Figure  5 
(V),  the  first  of  two  sawtooth  pulses  marks  the  correct  end  of  LVET. 

This  pulse  is  enabled  when  analog  processing  and  input  parameters  are 
selected  to  cause  the  program  to  identify  the  first  sawtooth  at  digital 
processing  time.  When  the  carotid  pulse  had  only  one  notch  and  a 
single  corresponding  sawtooth  pulse  was  present,  the  program  used  it 
as  LVET  termination.  A character  (#)  is  also  printed  on  the  printout's 
right  margin  to  indicate  which  channel  is  used  for  LVET  termination. 

B,  the  beginning  of  LVET,  is  located  in  digital  mode.  A search 
is  made  on  the  carotid  pulse  channe]  for  the  maximum  point,  X,  between 
two  consecutive  R waves  (Fig.  6).  The  minimum  point,  N,  is  located 
between  the  first  R wave  and  X.  Six  points  are  selected  between  N and 
X,  with  a least-squared  technic  used  to  find  line  1.  Line  2 is  a 
horizontal  line  intersecting  the  minimum  point,  N.  The  intersection 
of  these  two  lines  at  B is  then  computed  to  mark  the  beginning  of  LVET. 

The  STI  are  calculated  from  the  ST  from  single  cardiac  cycles. 

This  includes  Q-S2,  LVET,  Q-Si,  PEP  [(Q-S2)  - LVET],  and  IVCT  [(S^S^  - 
LVET],  The  STI  indices  are  computed  beat  by  beat,  using  the  standard 
Weissler  equations  shown  in  Table  1.  All  STI  and  STI  indices  are 
printed  in  milliseconds.  Measurements  of  the  beat  number,  the  heart  rate 
(beats/mln),  the  period  of  R-R  interval,  and  the  ratio  of  PEP  to  LVET 
are  also  computed  and  printed  as  shown  in  Figure  7.  At  the  top  of  each 
page,  the  program  prints  a heading  with  subject  name,  date  of  experiment, 
type  of  stress,  analog  tape  number,  and  time  code  (analysis  period). 

If  the  STI  had  a missing  component,  Sj  not  found,  asterisks  are  printed 
for  this  beat  in  the  column  requiring  Sj  (Q-Si,  S1-S2,  and  IVCT)  (Fig. 
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Figure  6.  The  R wave  sawtooth  pulses  are  used  to  bracket  the  events  in 


a cardiac  cycle.  Six  points  are  selected  between  N (minimum 
value)  and  X (maximum  value)  with  a least-squared  technic 
used  to  find  line  1.  The  intersection  of  line  1 and  line  2 
(horizontal  line  drawn  through  N)  marks  the  beginning  of 
LVET. 
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RESULTS 


The  STI  for  baseline  stress  conditions  were  manually  calculated  and 
compared  to  computerized  values.  In  all  baseline  and  stress  recovery 
cases,  the  automated  computations  (c.f.  Fig.  7)  agreed  with  manually 
reduced  values  within  2%  for  mean  levels.  Because  the  phonocardio- 
grams  and  carotid  pulse  signals  became  noisy  with  large  D.C.  artifacts 
during  stress  periods,  both  automated  and  manual  STI  were  inconsistently 
detected  during  stress  periods.  Nevertheless,  the  manual  and  automatic 
values  for  periods  of  reliable  signals  agreed  within  10%  in  all  cases. 

Because  the  beat-to-beat  variations  in  STI  stand  out  as  a unique 
advantage  of  the  automated  methodology,  mean  values  will  not  be  con- 
sidered here.  Rather,  the  instantaneous  values  will  be  presented  as 
depicted  in  Figure  7.  Figure  7 lists  reduced  STI  data  for  the  recovery 
period  immediately  following  +3  Gz  acceleration  exposure.  The  cardiac 
period  oscillates  between  values  of  approximately  750  and  950  msec, 
with  relatively  relaxed  breathing.  Computed  STI  are  listed  on  an 
individual  beat  basis  in  both  uncorrected  (PEP,  LVET,  Q-S2)  and  rate- 
corrected  (PEPI,  ETI,  Q-S2I)  form.  In  all  cases,  the  rate-corrected 
figures  are  larger  than  uncorrected,  as  expected  by  the  use  of  Weissler 
correction  factors.  Figures  8-10  display  a typical  set  of  baseline 
determinations  including  cardiac  period,  uncorrected  electromechanical 
systole,  and  heart  rate-corrected  preejection  period  as  plotted  on  a 
beat-to-beat  basis.  The  relatively  regular  oscillations  in  cardiac 
period  (Fig.  8)  follow  the  frequency  of  respiration,  apparently 
reflecting  normal  sinus  arrhythmia  and  central  blood  volume  changes 
during  normal  breathing.  Total  electromechanical  systole  (Fig.  9) 
closely  follows  the  same  oscillatory  pattern,  with  maximal  lengths  of 
systole  corresponding  to,  or  slightly  preceding,  maxima  for  the  cardiac 
period.  The  preejection  period  of  Fig.  10,  although  corrected  for 
changes  in  the  instantaneous  heart  rate,  similarly  varies  with  respira- 
tion. The  occurrence  of  maximal  preejection  periods,  however,  corre- 
sponds in  most  cases  to  minimal  cardiac  period  and  hence  to  reduced 
duration  of  total  electromechanical  systole.  At  beat  numbers  12,  20, 
and  26,  for  example,  the  PEP  is  at  high  peak  values  while  cardiac  period 
and  Q-S2  are  at  minima.  The  oscillatory  phenomena,  occurring  in 
response  to  normal  resting  respiration,  characterize  the  individual's 
cardiovascular  responses  to  the  rather  well-defined  mechanical  and 
neurologic  events  accompanying  lung  inspiration  and  expiration. 

In  Figure  1 1 similar  beat-to-beat  phenomena  are  demonstrated  for  a 
human  subject  recovering  from  acceleration.  Because  of  the  relatively 
high  heart  rate  and  respiration  rate  during  this  recovery  period,  the 
respiration- induced  oscillations  in  both  cardiac  period  (Fig.  11)  and 
LVET  (Fig.  12)  are  faster  with  relatively  sharper  peaks.  As  noted  for 
Q-S2  during  baseline  the  LVET  maxima  and  minima  are  in  phase  with  those 
for  cardiac  period.  The  data  plotted  in  Figures  8 through  12  are 
typical  for  all  subjects  measured  in  the  respect  that  respiration-induced 
oscillations  in  STI  intervals,  with  corresponding  relationships  between 
individual  intervals,  were  of  the  same  nature. 
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Figure  8.  Beat-by-beat  prestress  baseline  intervals 
for  human  subject  during  cardiac  period. 
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Figure  9.  Beat-by-beat  prestress  baseline  intervals  for  human  subject 
with  uncorrected  electromechanical  systole. 
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Figure  10.  Beat-by-beat  prestress  baseline  Intervals  for  human  subject 
with  electromechanical  systole  corrected  for  heart  rate. 
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DISCUSSION 


The  recording  and  analysis  of  systolic  time  intervals  provide 
clinicians  and  physiologic  researchers  with  a noninvasive,  atraumatic 
means  of  characterizing  the  cardiovascular  system's  condition.  The 
extensive  time  and  effort  required  for  manual  reduction  of  data  have 
precluded  the  routine  use  of  STI  analysis  in  clinical  or  experimental 
environments.  The  automated  methodology  detailed  in  this  report  has 
been  applied  to  human  subjects  during  various  combinations  of  environ- 
mental stress. 

STI  analysis  during  periods  of  stress  was  restricted  to  episodic 
information  because  of  mot ion- induced  noise  in  one  or  more  of  the  raw- 
signal  channels,  thereby  precluding  accurate  detection  of  some  of  the 
systolic  events  (see  noisy  pulse  in  Fig.  4).  Although  significant 
periods  of  transient  stress  response  are  thereby  lost,  sufficient  details 
of  the  integrated  physiologic  response  are  provided  during  periods 
Immediately  following  stress  exposure,  whenever  steady  state  conditions 
continue  to  prevail.  The  hybrid  system  has  further  demonstrated 
capabilities  not  found  in  other  existing  computer  programs:  (1) 
detections  of  individual  systolic  events  are  generally  mutually 
independent  because  of  the  separate  characterization  of  individual 
waveform  channels  during  analog  pattern  recognition;  (2)  the  analog 
computer  segment  provides  detection  capability  at  integral  multiples 
of  real  time,  thereby  greatly  reducing  time  for  processing;  (3)  long 
periods  of  time  are  reliably  analyzed  on  a beat-by-beat  basis  as  opposed 
to  an  averaged  scheme:  data  compression  via  sawtooth  technic  enables  the 
digital  program  to  process  fewer  discrete  samples  with  an  associated 
decreased  cost  for  input-output  and  central  processing  time. 

As  described  in  the  methodology  section,  certain  parameters  are 
input  into  the  digital  section  of  computer  analysis  on  an  individual 
subject  basis.  This  provision  assures  accurate  identification  of 
complicated  or  special  carotid  pulse  waveforms  and  thereby  takes  into 
consideration  individual  variations  among  subjects.  Only  minor  modlfica 
tions  are  required  for  conversion  of  the  technic  to  real-time  capability 
but  the  authors  feel  that  the  nature  of  phonocardiogram  and  carotid 
pulse  recordings  requires  special  analytical  preparation  and,  at  least 
in  the  research  setting,  large  quantities  of  information  could  be  lost 
or  deteriorated  by  inappropriate  analog  setting  occurring  in  real-time 
operation. 

Sample  data,  depicted  in  Figures  7 through  12,  demonstrate  beat-to- 
beat  variations  in  STI  that  oscillate  with  the  frequency  of  respiration. 
The  instantaneous  changes  in  preejection  period  (Fig.  10)  during  resting 
ventilation  reflect  the  neurological  and  mechanical  influence  of  lung 
inflation  and  deflation  on  left  heart  function.  In  particular,  the 
amplitude  of  PEP  oscillatory  change  represents  a quantitative  measure 
of  left  ventricular  filling  and  sinoatrial  node  electrical  activity. 

The  LVET  values,  during  recovery  from  acceleration  (Fig.  12),  further 
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demonstrate  two  simultaneous  dynamic  processes:  (1)  gradual  Increase 
in  ejection  time  over  the  complete  period  of  40  beats  corresponding  to 
an  increasing  cardiac  period  (Fig.  11),  with  concomitant  return  of  venous 
blood  to  the  central  volume;  (2)  oscillatory  changes  corresponding  to 
lung  inspiration  and  expiration.  These  early  differentiated  long-  and 
short-term  changes  in  LVET  reflect  the  underlying  cardiovascular  pro- 
cesses of  increasing  venous  return,  decreasing  systemic  resistance, 
increasing  central  blood  volume,  decreasing  myocardial  inotropy; 
and  reduced  ventilation.  The  relative  contributions  of  these  individual 
components  in  the  integrated  physiologic  response  can  be  quantitatively 
assessed  by  comparing  the  dynamic  changes  in  STI  among  one  another  and 
serially  evaluating  individual  values  with  reference  to  baseline  control. 

The  cardiovascular  effects  of  certain  pharmacological  agents  with 
positive  inotropic  effects,  for  instance,  could  be  evaluated  via  this 
automated  scheme  with  respect  to  individual  subject  response  by  assessing 
serial  changes  in  PEP,  LVET,  and  PEP/LVET.  While  the  time  and  effort 
required  for  manual  reduction  of  the  STI  in  such  a clinical  or  physiologic 
investigation  would  be  prohibitive,  the  application  of  the  automated 
method  converts  the  study  to  a routine  level.  The  method's  outstanding 
capability  to  process  long  periods  of  data  enables  a complete  cardio- 
vascular response  profile,  including  transients  and  steady  states. 

The  capability  of  the  automated  STI  analysis  to  provide  beat-to- 
beat  values  over  long  periods  of  time  extends  the  applications  and  values 
of  noninvasive  cardiac  records.  The  literature  to  date  has  reported  mean 
values  for  STI  (c.f.  ref.  5 and  15)  with  little  or  no  reference  to 
instantaneous  changes  and  their  significance.  The  method  reported  here, 
however,  has  evidenced  the  potential  informational  value  of  individual 
heart  beat  analysis.  Respiration  and  its  interaction  with  cardiac  and 
vascular  dynamics  has  been  shown  here  to  be  a useful  perturbation  when 
assessing  cardiac  status  via  STI.  Both  long-term  (multibeat)  and  short- 
term (beat-to-beat)  changes  in  systolic  intervals  can  be  accurately 
monitored  and  assessed  as  reflectors  of  cardiac  response  components 
(with  quantitation  ensured  by  the  computerized  scheme).  Investigators 
can  thus  ascertain  reliable  physiologic  and  pathologic  response  data 
with  noninvasive  means.  Minor  programming  changes  are  easily  accommodated 
with  this  methodology,  thus  ensuring  satisfaction  of  individual  investi- 
gator preference  for  special  analog  filtering  and  digital  computation 
correction. 
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